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Introduction 

The  plasma  population  in  the  inner  magnetosphere  has  two  sources,  namely  the 
dense,  cold  ionosphere  and  the  cool  solar  wind  [  Williams ,  1987].  The  relative  importance 
of  these  sources,  which  varies  as  a  function  of  geomagnetic  conditions,  can  be  identified 
from  higher-mass  ions:  Whereas  the  solar  wind  contains  alpha-particles  and  small 
amounts  of  heavier  elements,  the  ionosphere  provides  mostly  singly  charged  helium 
and  oxygen  in  addition  to  hydrogen.  The  relative  importance  of  these  sources  is  also 
a  function  of  location:  the  geopause  separates  the  domain  in  the  inner  magnetosphere 
dominated  by  ionospheric  source  from  the  solar  wind-dominated  outer  magnetosphere 
[Winglee,  1998]. 

The  outflow  of  ionospheric  ions  is  largely  controlled  by  the  thermospheric 
temperature,  which  varies  both  on  short  time  scales  as  a  response  to  geomagnetic 
activity  as  well  as  over  the  solar  cycle.  Both  magnetospheric  substorms  and  geomagnetic 
storms  heat  the  thermosphere  and  hence  lead  to  enhancement  of  the  ionospheric  ion 
outflow  from  the  dayside  cleft  ion  fountain  [ Cladis ,  1986]  and  from  the  nightside  auroral 
zone  [Daglis  et  al.,  1992].  On  the  other  hand,  increase  in  the  solar  UV  irradiance  can 
increase  the  thermospheric  temperature  by  a  factor  of  four  at  solar  maximum  relative 
to  the  temperature  during  solar  minimum  [Moore  et  al.,  1999],  which  leads  to  generally 
more  ionospheric  plasma-rich  magnetosphere  during  the  solar  maximum  [Young  et  al., 
1982}. 

The  inner  magnetosphere  consists  of  several  distinct  plasma  populations:  The 
cold  plasmasphere,  the  energetic  particles  in  the  radiation  belts,  the  hot  ring  current, 
and  the  inner  edge  of  the  plasma  sheet  are  all  partially  colocated.  The  ring  current  is 
especially  important  in  terms  of  magnetospheric  dynamics,  as  it  largely  controls  the 
inner  magnetospheric  magnetic  field  configuration.  Therefore,  it  is  especially  interesting 
to  examine  how  the  inner  magnetosphere  dynamics  depends  on  the  composition  of  the 
ring  current.  As  the  ring  current  ions  have  characteristic  energies  of  several  tens  of 


1 


Introduction 


The  plasma  population  in  the  inner  magnetosphere  has  two  sources,  namely  the 
dense,  cold  ionosphere  and  the  cool  solar  wind  [Williams,  1987].  The  relative  importance 
of  these  sources,  which  varies  as  a  function  of  geomagnetic  conditions,  can  be  identified 
from  higher-mass  ions:  Whereas  the  solar  wind  contains  alpha-particles  and  small 
amounts  of  heavier  elements,  the  ionosphere  provides  mostly  singly  charged  helium 
and  oxygen  in  addition  to  hydrogen.  The  relative  importance  of  these  sources  is  also 
a  function  of  location:  the  geopause  separates  the  domain  in  the  inner  magnetosphere 
dominated  by  ionospheric  source  from  the  solar  wind-dominated  outer  magnetosphere 
[Winglee,  1998]. 

The  outflow  of  ionospheric  ions  is  largely  controlled  by  the  thermospheric 
temperature,  which  varies  both  on  short  time  scales  as  a  response  to  geomagnetic 
activity  as  well  as  over  the  solar  cycle.  Both  magnetospheric  substorms  and  geomagnetic 
storms  heat  the  thermosphere  and  hence  lead  to  enhancement  of  the  ionospheric  ion 
outflow  from  the  dayside  cleft  ion  fountain  [ Cladis ,  1986]  and  from  the  nightside  auroral 
zone  [Daglis  et  ai,  1992],  On  the  other  hand,  increase  in  the  solar  UV  irradiance  can 
increase  the  thermospheric  temperature  by  a  factor  of  four  at  solar  maximum  relative 
to  the  temperature  during  solar  minimum  [ Moore  et  ai,  1999],  which  leads  to  generally 
more  ionospheric  plasma-rich  magnetosphere  during  the  solar  maximum  [Young  et  ai, 
1982], 

The  inner  magnetosphere  consists  of  several  distinct  plasma  populations:  The 
cold  plasmasphere,  the  energetic  particles  in  the  radiation  belts,  the  hot  ring  current, 
and  the  inner  edge  of  the  plasma  sheet  are  all  partially  colocated.  The  ring  current  is 
especially  important  in  terms  of  magnetospheric  dynamics,  as  it  largely  controls  the 
inner  magnetospheric  magnetic  field  configuration.  Therefore,  it  is  especially  interesting 
to  examine  how  the  inner  magnetosphere  dynamics  depends  on  the  composition  of  the 
ring  current.  As  the  ring  current  ions  have  characteristic  energies  of  several  tens  of 


1 


kilovolts,  active  processes  are  required  to  accelerate  the  ions  from  the  source  temperature 
to  the  observed  ring  current  temperatures  [ Cladis ,  1986,  Moore  et  a/.,  1999]. 

Geomagnetic  storms  are  created  by  a  variety  of  strong  disturbances  originating  from 
the  Sun.  For  example,  coronal  mass  ejections  (CME)  and  the  subsequent  interplanetary 
magnetic  clouds  are  characterized  by  strong  and  slowly  varying  interplanetary  magnetic 
field  (IMF).  When  the  IMF  points  southward,  CMEs  are  effective  drivers  of  strong 
geomagnetic  activity  [ Burlaga ,  1991].  One  of  the  main  characteristics  of  geomagnetic 
storms  is  the  formation  of  an  intense  ring  current  [ Hamilton  et  al .,  1988],  often  with 
an  enhanced  ionospheric  component  present  [Kremser  et  al ,  1993;  Roeder  et  a/.,  1999]. 
Although  strong  substorms  occur  during  geomagnetic  storms,  storms  are  more  than  a 
composite  of  several  substorms  [ McPherron ,  1997]. 

Both  geomagnetic  storms  and  substorms  are  rich  phenomena  associated  with  a 
number  of  different  processes  in  various  magnetospheric  regions.  In  order  statistically 
to  characterize  substorms  and  storms,  indices  derived  from  ground  magnetic  records  are 
often  used.  Substorms  are  mostly  effective  at  high  latitudes,  and  hence  the  substoim 
intensity  is  characterized  by  the  auroral  electrojet  indices  AU/AL  derived  from  auroral 
zone  magnetograms  or  the  3-hour  Kp  index  giving  a  more  time- averaged  view  of  the 
activity.  On  the  other  hand,  geomagnetic  storms  create  a  strong  ring  current  whose 
effects  can  be  measured  by  mid-latitude  magnetograms.  The  average  disturbance  o\er 
several  stations  is  called  the  Dst  index  and  is  used  to  characterize  the  storm  intensity. 

This  index  is  usually  given  as  an  hourly  average. 

The  storm  and  substorm  activity  vary  also  on  a  longer  time  scale  over  the  11-year 
solar  cycle.  Whereas  the  substorm  activity  often  shows  a  two-peaked  structure,  one 
slightly  before  the  solar  maximum  and  the  other  during  the  early  declining  phase 
[Nevanlinna  and  Pulkkinen ,  1998],  the  storm  occurrence  peaks  during  the  declining 
phase,  when  strong  geomagnetic  activity  is  driven  byr  coronal  mass  ejections  [Gosling 
et  al ,  1991]  and  high-speed  solar  wind  streams  originating  from  low-latitude  coronal 


2 


kilovolts,  active  processes  are  required  to  accelerate  the  ions  from  the  source  temperature 
to  the  observed  ring  current  temperatures  [Cladis,  1986,  Moore  et  al.,  1999], 

Geomagnetic  storms  are  created  by  a  variety  of  strong  disturbances  originating  from 
the  Sun.  For  example,  coronal  mass  ejections  (CME)  and  the  subsequent  interplanetary 
magnetic  clouds  are  characterized  by  strong  and  slowly  varying  interplanetary  magnetic 
field  (IMF).  When  the  IMF  points  southward,  CMEs  are  effective  drivers  of  strong 
geomagnetic  activity  [Burlaga,  1991].  One  of  the  main  characteristics  of  geomagnetic 
storms  is  the  formation  of  an  intense  ring  current  [Hamilton  et  al. ,  1988],  often  with 
an  enhanced  ionospheric  component  present  [Kremser  et  al.,  1993;  Roeder  et  al.,  1999]. 
Although  strong  substorms  occur  during  geomagnetic  storms,  storms  are  more  than  a 
composite  of  several  substorms  [McPherron,  1997]. 

Both  geomagnetic  storms  and  substorms  are  rich  phenomena  associated  with  a 
number  of  different  processes  in  various  magnetospheric  regions.  In  order  statistically 
to  characterize  substorms  and  storms,  indices  derived  from  ground  magnetic  records  are 
often  used.  Substorms  are  mostly  effective  at  high  latitudes,  and  hence  the  substorm 
intensity  is  characterized  by  the  auroral  electrojet  indices  AU/AL  derived  from  auroral 
zone  magnetograms  or  the  3-hour  Kp  index  giving  a  more  time-averaged  view  of  the 
activity.  On  the  other  hand,  geomagnetic  storms  create  a  strong  ring  current  whose 
effects  can  be  measured  by  mid-latitude  magnetograms.  The  average  disturbance  over 
several  stations  is  called  the  Dst.  index  and  is  used  to  characterize  the  storm  intensitv. 
This  index  is  usually  given  as  an  hourly  average. 

The  storm  and  substorm  activity  vary  also  on  a  longer  time  scale  over  the  11-year 
solar  cycle.  Whereas  the  substorm  activity  often  shows  a  two-peaked  structure,  one 
slightly  before  the  solar  maximum  and  the  other  during  the  early  declining  phase 
[Nevanlinna  and  Pulkkinen ,  1998],  the  storm  occurrence  peaks  during  the  declining 
phase,  when  strong  geomagnetic  activity  is  driven  by  coronal  mass  ejections  [Gosling 
et  al.,  1991]  and  high-speed  solar  wind  streams  originating  from  low-latitude  coronal 


2 


holes  [Feldman  et  al.,  1978].  The  11-year  cycle  also  causes  other  variations  in  the  solar 
effects  on  Earth,  one  of  the  most  important  for  the  ionospheric  processes  being  the 
increase  of  the  UV  flux  with  solar  activity,  which  is  a  potentially  important  part  of  the 
ionosphere-magnetosphere  coupling.  The  F10.7  radio  flux  is  often  used  as  a  proxy  for 
the  UV  flux  variations. 

The  first  three  years  of  operation  of  the  POLAR  spacecraft  allow  for  detailed 
examination  of  the  plasma  composition  during  a  very  quiet  period  near  the  solar 
minimum  and  during  the  rising  phase  of  the  solar  cycle.  The  Charge  and  Mass 
Magnetospheric  Ion  Composition  Measurement  /  Magnetospheric  Ion  Composition 
Sensor  (CAMMICE/MICS)  [  Wilken  et  al. ,  1992]  is  capable  of  measuring  ions  roughly 
in  the  energy  range  from  1  keV  to  200  keV,  and  separating  particle  species  as  well  as 
their  charge  states.  The  highly  elliptical  polar  orbit  with  an  18-hour  orbital  period 
compromises  time  resolution,  but  allows  for  observations  at  all  magnetic  L  shells  above 

L  =  3. 

In  this  paper,  we  compare  the  ion  flux  measurements  with  the  interplanetary 
magnetic  field  and  solar  wind  measurements  and  with  the  geomagnetic  activity  indices. 
Both  average  flux-values  and  average  energy  density  values  are  computed  for  hydrogen, 
singly  and  doubly  charged  helium,  and  for  singly  charged  oxygen.  The  IMF,  solar 
wind,  and  activity  parameters  are  averaged  over  the  periods  when  the  measurements 
were  made,  and  the  correlation  of  the  flux  and  energy  density  values  with  both  the 
interplanetary  drivers  and  the  proxies  for  geomagnetic  activity  are  examined.  A  specific 
goal  of  the  study  is  to  examine  the  ion  composition  changes  during  geomagnetically 
active  periods  (especially  storms)  during  different  phases  of  solar  activity. 

Of  the  ion  species  selected  for  examination,  H+  originates  both  from  the  ionosphere 
and  the  solar  wind,  whereas  0+  comes  only  from  the  ionosphere  and  He+~  only  from 
the  solar  wind  On  the  other  hand,  He+  can  be  either  of  ionospheric  or  of  solar  wind 
origin  (He++  charge-exchanged).  Therefore,  enhancements  of  He+  and  He*-  need  to  be 
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in  the  energy  range  from  1  keV  to  200  keV,  and  separating  particle  species  as  well  as 
their  charge  states.  The  highly  elliptical  polar  orbit  with  an  18-hour  orbital  period 
compromises  time  resolution,  but  allows  for  observations  at  all  magnetic  L  shells  above 
L  =  3. 

In  this  paper,  we  compare  the  ion  flux  measurements  with  the  interplanetary 
magnetic  field  and  solar  wind  measurements  and  with  the  geomagnetic  activity  indices. 
Both  average  flux-values  and  average  energy  density  values  are  computed  for  hydrogen, 
singly  and  doubly  charged  helium,  and  for  singly  charged  oxygen.  The  IMF,  solar 
wind,  and  activity  parameters  are  averaged  over  the  periods  when  the  measurements 
were  made,  and  the  correlation  of  the  flux  and  energy  density  values  with  both  the 
interplanetary  drivers  and  the  proxies  for  geomagnetic  activity  are  examined.  A  specific 
goal  of  the  study  is  to  examine  the  ion  composition  changes  during  geomagnetically 
active  periods  (especially  storms)  during  different  phases  of  solar  activity. 

Of  the  ion  species  selected  for  examination,  H+  originates  both  from  the  ionosphere 
and  the  solar  wind,  whereas  0+  comes  only  from  the  ionosphere  and  He+~  only  from 
the  solar  wind.  On  the  other  hand,  He+  can  be  either  of  ionospheric  or  of  solar  wind 
origin  (He++  charge-exchanged).  Therefore,  enhancements  of  He+  and  He~~  need  to  be 
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examined  together;  here  we  do  not  attempt  to  distinguish  the  portion  of  He+  coming 
the  charge-exchange  processes  [Kremser  et  ai,  1993]. 


Data  set 

The  data  used  in  this  study  comprises  pitch-angle-averaged  fluxes  measured  by  the 
CAMMICE/MICS  instrument  [Wilken  et  al.,  1992]  onboard  POLAR  during  September 
1,  1996,  -  March  31,  1999.  MICS  measures  a  full  distribution  function  in  96  seconds, 
which  is  why  the  data  were  averaged  to  that  resolution.  The  data  are  tagged  with 
POLAR  orbital  information  and  L  shell  values  obtained  from  the  POLAR  data  base. 
Note  that  the  L  values  do  not  account  for  the  actual  state  of  the  magnetosphere, 
which  can  introduce  errors  in  the  L  value  estimations  during  highly  disturbed  periods. 
However,  accounting  for  the  actual  mapping  was  not  possible,  because  this  hass  to  be 
done  on  an  event-by-event  basis,  and  because  there  are  no  good  magnetic  field  models 
for  the  storm-time  inner  magnetosphere. 

The  measured  flux,  j(E,L),  was  integrated  over  all  measured  energies 

rErnax 

J(L)=  dEj(E,L)  (1) 

''Emin 

to  give  the  total  flux,  J(L),  as  a  function  of  L.  Furthermore,  we  compute  the  average 
and  peak  fluxes  between  L  =  3  and  L  =  8  for  each  pass.  In  the  computations,  the 
particle  pitch-angle  distributions  were  assumed  to  be  sufficiently  isotropic  that,  within 
this  L  range,  no  corrections  to  equatorial  flux  values  were  necessary.  Figure  la  shows  a 
sample  of  the  flux  values  as  a  function  of  L  on  Oct  23,  1996.  Both  the  maximum  and 
average  fluxes  are  marked  in  the  figure  for  hydrogen,  helium,  and  oxygen. 

If  a  pitch-angle  correction  were  to  be  applied,  the  standard  way  would  be  to  assume 
a  dependence  jeq(E,a )  =  jmeas(E,  n/2)  sin71  a,  where  jeq  and  jmeas  are  the  equatorial 
and  measured  fluxes,  and  a  is  the  pitch-angle.  However,  detailed  examinations  of  the 
distribution  functions  suggested  that,  for  the  most  part,  the  distributions  within  the 
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Figure  la.  POLAR/CAMMICE/MICS  observations  of  H  ,  He  ,  He  ,  and  O  ions  during  a  geomag¬ 
netic  storm  on  Oct  23,  1996.  The  left  panel  shows  the  flux  values,  whereas  the  right  panel 
shows  the  energy  densities  as  a  function  of  L  shell. 


fyjjQg  energy  range  were  Quite  isotropic,  and  that  a  common  anisotropy  factor  (rz)  could 
not  be  found  for  the  entire  data  set.  Because  finding  the  correct  n  is  quite  a  complicated 
process  [Roeder  et  al. ,  1999],  we  assumed  n  =  0  in  the  statistical  examination.  The 
result  of  this  approximation  leads  to  a  slight  underestimation  of  the  fluxes  (especially 
for  heavier  ions)  during  storm  periods  and/or  during  lower  solar  activity  when  the 
anisotropy  is  likely  to  be  larger.  On  the  other  hand,  n  typically  becomes  large  only  for 
energies  much  higher  than  those  covered  by  the  CAMMICE/MICS  instrument. 

Figure  lb  shows  an  overview  of  the  period  in  question.  The  top  panels  show  four 
channels  of  the  MICS  instrument,  RO  (H+),  Rl  (He+),  R2  (He++),  and  R3  (0+)  in  an 
L  shell  -  time  format,  where  the  total  flux  intensity  is  color  coded.  The  following  two 
panels  show  two  flux  ratios,  0+/H+  and  0+/He+.  The  bottom  panels  show  the  Dst 
index  characterizing  the  geomagnetic  activity  conditions,  and  the  FI  0.7  flux  indicating 
the  level  of  solar  activity.  Storm  periods  (Dst  <  —50  nT)  are  marked  red  in  the 
Dst  panel.  The  F10.7  flux  clearly  shows  how  the  solar  activity  starts  to  rise  around 
September  1997  (Day  600  in  Figure  lb).  A  similar  trend  of  rising  activity  is  also  evident 
in  the  Dst  data,  with  more  intense  storm  activity  after  September  1997  than  during  the 
solar  minimum  period  (September  1996  -  September  1997). 

The  particle  fluxes  show  the  bulk  of  the  ring  current  particles  mainly  in  the  region 
between  L  =  4  and  L  -  6.  All  species  show'  periodic  enhancements,  which  are  coincident 
with  the  Dst  enhancements  and  hence  associated  w'ith  geomagnetic  storms.  Both  the 
solar  wind  origin  He++  and  ionospheric  0+  are  enhanced  during  storms,  although  the 
effects  in  He++  are  weaker  and  the  population  is  generally  at  higher  L  values.  Note 
also  how  the  general  level  of  oxygen  increases  from  the  solar  minimum  toward  higher 
solar  activity.  The  panel  showing  the  0+/H+  ratio  show's  that  the  oxygen  flux  is  low- 
relative  to  the  hydrogen  flux,  typically  only  a  few  percent.  How-ever,  the  amount  of 
oxygen  increases  with  the  solar  activity,  but  even  at  the  end  ofthe  period,  the  ratio 
rarely  reaches  10%.  Similarly,  the  Cr/He4'  show's  that  the  helium  fluxes  dominate  over 
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in  the  Dst  data,  with  more  intense  storm  activity  after  September  1997  than  during  the 
solar  minimum  period  (September  1996  -  September  1997). 

The  particle  fluxes  show  the  bulk  of  the  ring  current  particles  mainly  in  the  region 
between  L  =  4  and  L  =  6.  All  species  show  periodic  enhancements,  which  are  coincident 
with  the  Dst  enhancements  and  hence  associated  with  geomagnetic  storms.  Both  the 
solar  wind  origin  He++  and  ionospheric  0+  are  enhanced  during  storms,  although  the 
effects  in  He++  are  weaker  and  the  population  is  generally  at  higher  L  values.  Note 
also  how  the  general  level  of  oxygen  increases  from  the  solar  minimum  toward  higher 
solar  activity.  The  panel  showing  the  0+/H+  ratio  shows  that  the  oxygen  flux  is  low 
relative  to  the  hydrogen  flux,  typically  only  a  few  percent.  However,  the  amount  of 
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Figure  lb.  POLAR/CAMMICE/MICS  measurements  of  the  integrated  flux  vs  time  and  L-shell  over 
the  period  from  September  1,  1996  to  March  31,  1999.  Panels  from  top  to  bottom:  H+, 
Hc+,  He++,  0+,  flux  ratios  0+/H+,  and  0+/Hc+.  Note  that  the  scale  for  H+  is  different 
from  the  other  ion  species.  The  two  lowest  panels  show  the  Dst  index  and  the  FI 0.7  flux. 
The  time  is  given  in  days  of  the  year  1996. 


the  oxygen  fluxes,  which  are  typically  less  than  50%  of  the  helium  fluxes.  However,  as 
the  solar  activity  increases  during  1998  and  1999,  the  oxygen  fluxes  can  at  times  rise 
even  over  the  helium  flux  values. 

For  the  statistical  examination,  the  mean  total  flux  ( Jave )  (within  the  L  shell  range 
between  L  =  3  and  L  =  8),  the  peak  total  flux  (Jma *),  and  the  L  shell  where  the  peak 
flux  was  detected  ( Ljtmax )  were  tabulated.  The  L  value  limits  were  chosen  to  restrict  the 
examination  only  to  the  ring  current  region  -  high  energetic  ion  fluxes  can  occasionally 
be  found  near  the  dayside  magnetopause,  but  those  are  not  included  in  this  study. 

The  energy  density  per  unit  volume  associated  with  the  ring  current  ions  was 
computed  from 

w(L)  -  27Tyj2mq  dEVEj(E,L)  (2) 

^  Emin 

where  w  is  the  energy  density,  m  the  ion  mass,  q  the  charge  state,  and  the  integral  is 
over  energy  per  unit  charge.  Also  in  this  case,  the  mean  energy  density  (wave),  peak 
energy  density  (wmax),  and  the  L  shell  where  the  peak  energy  density  was  detected 
(LWimaX)  were  tabulated. 

All  these  observations  were  tagged  with  geomagnetic  activity  parameters,  solar 
wind,  and  interplanetary  magnetic  field  information.  The  parameters  include  the 
planetary  activity  index  Kp,  the  Dst  index,  solar  wind  velocity  (V),  and  the  north-south 
component  of  the  interplanetary  magnetic  field  (Bz).  Furthermore,  to  examine  the 
solar  wind  -  magnetosphere  coupling,  we  computed  the  e  parameter  [Perreault  and 
Akasofu,  1978]  for  these  periods.  All  these  parameters  were  averaged  over  the  period 
during  which  POLAR  passed  through  the  ring  current  region  and  measured  the  flux 
and  energy  density  values.  Because  of  the  18-hour  orbit  of  POLAR,  it  passes  through 
the  ring  current  region  relatively  quickly,  which  suggests  that  for  the  most  part,  large 
temporal  changes  do  not  occur  during  any  one  pass.  On  the  other  hand,  due  to  the 
large  absences  of  data  at  low  L  values  when  POLAR  is  near  apogee,  it  is  impossible  to 
determine  the  exact  time  of  formation  of  the  measured  particle  populations.  Therefore. 
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any  time  delays  in  activity  parameters  have  not  been  included.  Furthermore,  12-hour 
averages  of  these  parameters  were  computed  such  that  the  12-hour  period  ended  at 
the  time  when  POLAR  completed  its  ring  current  crossing.  Thus,  these  average  values 
represent  the  parameters  prior  to  and  during  the  crossings.  Finally,  the  minimum  of  Dst 
and  the  maximum  of  Kp  during  each  of  the  12-hour  intervals  were  determined  and  used 
in  the  correlation  analysis. 

Statistical  results 

In  order  to  get  an  overview  of  the  variations  of  the  ion  composition  in  the 
inner  magnetosphere,  the  fluxes  and  energy  densities  were  evaluated  for  each  inner 
magnetosphere  crossing  during  the  time  period  from  September  1,  1996  to  March  31, 
1999.  The  bottom  panel  of  Figure  lb  shows  that  the  data  set  can  be  divided  into  three 
almost  equally  long  periods  characterized  by  “low  solar  activity”  (F10.7  <  100,  DOY 
245-600),  “rising  solar  activity”  (F10.7  ~  100,  DOY  600-950)  and  “high  solar  activity" 
(F10.7  >  100,  DOY  950-1186).  Differences  in  these  periods  thus  reflect  the  changes 
caused  by  the  longer-term  variations  in  the  solar  activity. 

Figure  2  shows  a  scatter  plot  where  the  hydrogen  energy  densities  have  been 
correlated  with  minimum  of  Dst,  average  Kp,  average  V ,  average  Bz ,  average  e,  and 
the  F10.7  flux.  The  Dst  minimum  and  all  averages  were  computed  during  the  12-hour 
averaging  interval  defined  above.  The  three  different  periods  have  been  marked  with 
different  symbols,  squares  for  DOY  245-600,  circles  for  DOY  600-950,  and  triangles 
for  DOY  950-1186.  Generally,  the  correlations  are  not  very  good,  the  best  ones  being 
the  minimum  of  Dst  during  the  prior  12  hours  and  the  average  of  Kp  during  the  prior 
12-hour  period.  The  different  years  do  not  show  much  variability  except  that  the 
highest  F10.7  values  occur  during  the  last  period,  which  was  how  the  data  set  was 
divided.  The  solar  wind  velocity  gave  somewhat  better  correlation  than  the  IMF  Bz. 
but  neither  correlation  is  above  R  ~  0.5.  Here  also  the  12-hour  average  values  give 
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Scatter  plots  of  H+  energy  density  (in  keV/cirf )  vs  the  minimum  of  Dst,  average  of  Kp, 
solar  wind  velocity,  IMF  Br  £  parameter,  and  the  FI 0.7  flux.  The  different  years  are 
shown  with  different  symbols  (squares  for  1996-1997,  circles  for  1997-1998,  and 
triangles  for  1998-1999). 


better  correlations  than  the  instantaneous  values  during  the  measurements,  but  the 
differences  are  not  large.  Surprisingly,  neither  the  e  parameter  measuring  the  efficiency 
of  the  solar  wind  entry  nor  the  F10.7  measuring  the  ionospheric  upflow  activity  show 
any  correlation  with  the  hydrogen  energy  density.  For  all  these  results,  the  correlations 
are  slightly  better  for  the  energy  density  values  than  for  the  averaged  flux  values. 

The  energy  densities  for  the  other  species  are  of  course  smaller  than  for  the 
dominant  hydrogen,  but  can  be  at  times  quite  significant,  as  demonstrated  in  Figure 
3.  This  shows  the  statistics  for  the  minimum  of  Dst,  the  average  of  Kp ,  the  average 
of  solar  wind  velocity,  and  F10.7  flux  for  He+,  He++,  and  0+.  For  the  Dst  index, 
the  correlation  is  best  for  oxygen,  indicating  its  importance  in  geomagnetic  storms 
which  will  be  discussed  in  more  detail  below.  Surprisingly,  the  best  correlation  with 
geomagnetic  activity  (Kp)  is  given  by  the  doubly  charged  (solar  wind  origin)  helium. 
This  is  probably  due  to  the  inward  penetration  of  the  plasma  sheet,  where  the  solar  wind 
helium  has  access  to.  This  might  be  taken  to  indicate  that  He++  would  correlate  well 
with  the  e  parameter,  but  this  is  not  the  case;  the  correlations  are  quite  low,  similar  to 
those  shown  for  H+  (not  shown).  On  the  other  hand,  the  solar  wind  velocity  correlation 
is  best  for  helium.  Note  also  that  the  range  of  energy  density  values  for  low  activity 
levels  is  much  wider  for  0+  than  for  the  other  species.  This  also  shows  a  clear  solar 
activity  dependence:  the  low  energy  density  levels  tend  to  concentrate  on  the  period  of 
low  solar  activity.  The  F10.7  flux  correlations  are  poor  in  all  cases,  but  are  somewhat 
larger  for  oxygen  (0.33)  than  for  the  other  species  (<  0.1),  which  is  one  illustration  of 
the  increase  of  the  oxygen  energy  densities  as  the  solar  activity  increases. 

It  is  interesting  to  note  that  all  the  correlations,  if  computed  for  the  three 
time  periods  separately,  show  best  correlations  for  the  period  of  rising  solar  activity 
(September  1997  -  September  1998).  The  reason  for  this  is  not  clear,  but  it  is  quite 
consistent  in  correlations  with  all  the  driver  and  activity  parameters,  both  for  the  energy 
densities  and  for  the  averaged  flux  values. 
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The  long-term  energy  density  content  of  the  ring  current  is  examined  in  Figure  4 
in  terms  of  27-day  running  averages  of  the  measured  fluxes.  The  top  four  panels  show 
each  of  the  species  separately,  whereas  the  ratios  of  the  heavier  elements  to  hydrogen 
and  to  singly  charged  helium  are  plotted  in  the  two  bottom  panels.  It  can  be  seen  that 
both  He+/H+  and  0+/H+  ratios  increase  as  the  solar  activity  increases,  but  that  the 
increase  in  the  oxygen  energy  density  is  much  more  significant.  The  0+  energy  density 
is  typically  around  3%  during  the  solar  minimum,  whereas  it  increases  to  around  10% 
during  higher  solar  activity.  This  is  evident  also  from  the  flux  values  in  the  fourth 
panel  as  well  as  from  the  bottom  panel  showing  the  increasing  significance  of  oxygen 
as  the  activity  incrases.  The  energy  densities  of  0+  and  He+  are  almost  equal  during 
the  solar  minimum  period,  but  for  the  higher  activity  period  the  0+  energy  density  is 
about  twice  the  He+  energy  density,  even  though  the  fluxes  are  smaller,  as  was  shown 
in  Figure  lb.  The  He++  to  H+  ratio  did  not  increase  substantially,  and  the  portion  of 
He++  remained  small  in  the  inner  magnetosphere  throughout  the  period.  The  vertical 
bars  show  the  standard  errors  (standard  deviation  divided  by  the  square  root  of  the 
number  of  samples,  shown  for  every  10th  data  point  to  avoid  crowding  of  the  plot). 

Storm  periods 

As  the  largest  enhancements  of  the  ionospheric  plasma  component  occur  during 
geomagnetic  storms,  the  entire  data  set  was  searched  to  examine  these  periods 
separately.  Storms  were  defined  as  periods  with  minimum  of  Dst  less  than  -50  nT,  and 
were  selected  by  manual  inspection  of  the  Dst  records  available  from  the  World  Data 
Center  in  Kyoto.  In  total,  44  storm  events  were  included  in  the  study.  Some  periods 
with  extended  Dst  activity  around  -50  nT  were  excluded  from  the  study,  because  they 
did  not  show  a  clear  peak  or  typical  storm-time  behavior.  Similarly,  peaks  during  the 
recovery  phases  of  storms  were  not  counted  as  separate  events. 

As  an  example,  Figure  5  shows  energy  -  L  spectrograms  from  several  consecutive 
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27-day  running  average  energy  densities 


Figure  4.  27-day  averaged  values  of  the  energy  densities  for  H+,  He+,  He+~\  and  0+.  The  two 

bottomlands  show  ratios  of  these  quantities  for  0+/H+,  He+/H+,  H++/H+,  and  0+/He+ 
and  He  /He  .  The  vertical  bars  show  the  standard  errors. 
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Figure  5.  Dst  index  during  the  storm  on  May  15,  1997.  POLAR/MICS  energy — L-shell  spectrograms 
during  a  storm  on  May  15,  1997  for  hydrogen  (top  six  panels)  and  for  oxygen  (bottom  six 
panels).  The  times  of  the  observation  periods  are  shown  with  red  bars  in  the  Dst  plot. 
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passes  during  the  May  15,  1997  storm  period.  The  top  panels  shows  the  Dst  index. 
The  red  lines  mark  the  periods  during  which  the  data  shown  in  the  panels  below  were 
recorded.  The  bottom  panels  show  the  spectrograms  in  a  format  where  the  dayside 
(1000  MLT)  passes  are  on  the  left  and  the  nightside  (2200  MLT)  passes  are  on  the 
right.  The  top  six  panels  show  the  hydrogen  fluxes  and  the  bottom  six  panels  show  the 
oxygen  fluxes.  Notice  the  first  strong  enhancement  at  low  energies  during  the  storm 
main  phase,  whereas  the  strongest  intensification  at  ring  current  energies  was  reached 
after  the  Dst  minimum  during  the  storm  recovery  phase. 

The  0+  behavior  shows  some  similarities  but  also  differences  to  the  H+  spectra. 
During  DOY  134  there  is  very  little  0+  in  the  magnetosphere.  Shortly  after  the  storm 
onset  (first  pass  on  DOY  135),  0+  appears  with  the  intensification  of  low-energy  H~, 
but  extending  to  much  higher  energies.  The  lack  of  0+  at  low  energies  is  probably  due 
to  the  low  sensitivity  of  the  instrument  in  this  energy  range.  Furthermore,  in  the  next 
dayside  pass,  it  seems  that  the  0+  population  is  located  closer  to  the  Earth  than  the 
lower-energy  H+  (tens  of  keV),  but  further  away  than  the  higher-energy  H+  (around 
and  above  100  keV),  which  extends  all  the  way  to  the  innermost  distance  covered  by 
the  measurements.  The  different  locations  of  the  H+  and  0+  populations  might  be 
indicative  of  separate  acceleration  processes  at  different  distances  or  different  injection 
locations  in  the  tail. 

In  the  statistical  examination,  we  consider  the  first  POLAR  crossings  through  the 
ring  current  region  after  the  minimum  of  Dst  was  reached,  i.e.,  measurements  that  were 
made  as  close  as  possible  to  the  ring  current  maximum  intensity.  In  the  event  shown 
above,  this  would  be  the  pass  well  into  the  storm  recovery.  This  outlines  one  of  the 
problems  in  this  correlation  study:  data  from  the  maximum  activity  periods  are  not 
always  available  in  this  data  set. 

Figure  6  shows  the  correlations  between  the  H+  energy  density  and  the  minimum 
of  Dst,  average  of  I<p.  average  of  the  solar  wind  speed,  and  the  F10.7  flux.  The  data  set 


16 


ave(V)  min(Dst) 


Wave(H+) 


Wave(H+) 


Wave(H+)  Wave(H+) 

Figure  6.  Scatter  plots  of  H+  energy  density  vs  the  minimum  of  Dst,  average  of  Kp,  solar  wind 

velocity,  and  the  FI 0.7  flux  for  44  storms  during  the  examined  time  interval.  The  different 
years  are  shown  with  different  symbols  (squares  for  1996-1997,  circles  for  1997-1998, 
and  triangles  for  1998-1999). 
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shows  only  the  measurements  closest  to  (but  after)  the  storm  maximum  as  defined  from 
the  minimum  of  Dst.  The  Dst  plot  still  shows  quite  large  variability,  and  the  correlation 
is  actually  slightly  poorer  than  it  is  for  the  full  dat  set.  On  the  other  hand,  correlation 
with  Kp  is  better  than  for  the  full  data  set.  It  is  also  clearly  seen  that  the  data  from 
the  high  solar  activity  period  (triangles)  make  the  correlation  worse;  these  periods  tend 
to  be  characterized  by  lower  fluxes  than  the  other  periods. 

The  picture  for  the  other  particle  species  becomes  much  more  organized  when  only 
the  storm  periods  are  considered.  Figure  7  shows  a  plot  similar  to  Figure  3,  where  the 
helium  and  oxygen  energy  densities  are  correlated  with  the  minimum  of  Dst,  average  of 
Kp,  and  average  of  solar  wind  speed.  The  correlations  are  clearly  higher  in  almost  all 
cases  than  for  the  full  data  set.  Especially,  the  0+  energy  density  shows  good  correlation 
with  the  Dst  minimum. 

Even  though  the  activity  parameters  show  a  clear  correlation  with  the  ring  current 
ion  fluxes,  the  lack  of  high  correlation  coefficients  points  out  that  the  processes  involve 
multiple  aspects,  and  that  the  resulting  ion  composition  is  a  function  of  several  factors. 
Figure  8  shows  one  more  parameter,  the  minimum  of  the  AL  index  during  the  storm 
period  correlated  with  the  average  fluxes  near  the  storm  maximum.  The  correlations 
are  not  very  good  for  H+  or  He,  but  it  is  comparable  to  the  other  activity  parameters 
for  0+. 

Figure  9  shows  a  superposed  epoch  analysis  of  the  44  storm  events.  The  top  panel 
shows  the  Dst  index,  the  three  following  panels  the  energy  densities  for  H+.  0+,  He4" 
and  He++.  Note  the  two-step  recovery  process  in  the  Dst-index,  where  the  more  rapid 
decay  after  the  main  phase  maximum  is  followed  by  a  slower  decay  lasting  over  the 
entire  8-day  period  shown  here.  The  bottom  panel  shows  the  L  values  where  the  peak 
energy  densities  were  recorded.  It  is  evident  from  this  plot  that  all  energy  density  levels 
rise  rapidly  after  the  main  phase  onset,  and  reach  maximum  values  by  the  Dst-minimum 
time.  Whereas  the  H+  energy  density  stays  high  throughout  the  8  days  shown  in  this 
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urc  7.  Scatter  plot  of  He+,  He++,  and  0+  energy  densities  vs  the  minimum  of  Dst,  average  of  Kp 
and  average  of  the  solar  wind  velocity  for  44  storm  events  during  the  examined  time 
interval.  Different  years  are  shown  with  different  symbols  (squares  for  1996-1997,  circles 
for  1997-1998,  and  triangles  for  1998-1999). 
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Figure  8.  Scatter  plot  of  He+,  He++,  and  0+  energy  densities  vs  the  minimum  of  AL  index  for  44 
storm  events  during  the  examined  time  interval.  Different  years  are  shown  with  different 
symbols  (squares  for  1996-1997,  circles  for  1997-1998,  and  triangles  for  1998-1999). 
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Figure  9. 


density  of  H  ,  He  ,  He  ,  and  O  ,  and  thcr  L-shcll  of  the  peak  flux  for  H  ,  He  ,  He 
and  O  .  The  pre-storm  levels  arc  indicated  by  the  horizontal  lines. 
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plot,  0+,  He+,  and  He++  decay  close  to  the  pre-storm  values  shown  by  the  horizontal 
solid  lines.  On  the  other  hand,  it  is  also  clear  that  even  helium  and  oxygen  decay  slowly 
enough  so  that  they  cannot  account  for  the  two-step  recovery  phase  where  the  first, 
rapid  recovery  lasts  for  less  than  one  day  and  the  slower  recovery  may  last  from  several 
days  to  two  or  three  weeks.  The  behavior  of  the  flux-values  are  virtually  identical  to  the 
energy  density  values  shown  in  Figure  9  (not  shown). 

The  vertical  bars  in  all  panels  show  the  standard  errors  in  the  averaging  procedure. 
For  the  energy  density  values,  the  standard  deviations  are  for  the  most  part  small 
enough  that  they  cannot  be  distinguished  from  the  plotting  symbols.  The  bottom  panel 
shows  the  general  inward  motion  of  the  ring  current  at  substorm  onset.  Whereas  the 
hydrogen  energy  density  maximum  stays  close  to  the  Earth  at  slightly  outside  L  =  4 
after  the  inward  injection,  the  maximum  energy  densities  of  oxygen  and  helium  return 
to  their  average  higher  L  values  within  the  same  time  frame  as  the  fluxes  decay. 

Comparison  of  the  Dst  profile  with  the  energy  density  variations  shows  that  H,  He-, 
and  0+  all  give  a  correlation  coefficient  of  about  0.90,  whereas  the  correlation  coefficient 
for  He++  is  somewhat  lower,  0.84  (not  shown).  Similar  comparison  of  the  L  values  of 
the  peak  energy  densities  with  the  Dst-profiles  gives  somewhat  lower  correlations,  for  H 
only  0.37,  for  0+  0.78,  for  He+  0.87,  and  for  He++  0.81. 

A  clear  distinction  in  the  correlation  coefficients  of  the  ionospheric  0+  and  solar 
wind  origin  He++  would  indicate  a  possible  cause  for  the  observed  two-step  recovery 
in  the  Dst-index:  The  L  value  changes  of  the  He++  energy  density  maximum  would 
then  imply  inward  penetration  and  consequent  retreat  to  higher  L  values  of  the  plasma 
sheet.  As  the  He  correlation  coefficients  are  somewhat  higher  than  those  for  0+.  this  is 
possibly  the  case,  but  further  case  studies  are  necessary  to  confirm  this  result. 
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Discussion 


We  have  shown  in  this  paper  that  the  ring  current  ion  composition  varies 
significantly  from  the  conditions  at  solar  minimum  (1996)  to  close  to  the  solar  maximum 
(1999):  The  average  0+  energy  density  increases  by  about  a  factor  of  5  during  the  rising 
phase  of  the  solar  cycle  from  the  minimum  values  in  1996,  and  the  more  frequent  storms 
under  these  conditions  give  rise  to  large  peaks  in  the  oxygen  and  helium  concentration 
in  the  inner  magnetosphere.  On  the  other  hand,  the  average  values  of  H+  and  He  energy 
density  show  variability  but  no  consistently  increasing  trend.  The  0+  flux  is  small 
(below  10%)  compared  to  the  hydrogen  flux,  and  the  average  energy  density  ranges 
from  a  few  %  at  solar  minimum  to  about  10%  at  high  solar  activity  time  in  early  1999. 
The  0+  flux  is  typically  smaller  than  the  Hc+  flux,  reaching  to  comparable  values  only 
during  the  latter  part  of  the  period  when  the  solar  activity  increased.  Analogously,  the 
energy  densities  of  0+  and  He+  are  about  equal  during  1996  and  1997,  whereas  the  O- 
energy  density  is  about  twice  the  He+  energy  density  during  the  higher  solar  activity 
period  in  1998  and  early  1999. 

The  correlation  with  the  geomagnetic  activity,  solar  wind  parameters,  and  overall 
solar  activity  is  very  complex.  Generally  the  ion  fluxes  are  poorly  correlated  with  the 
e  parameter,  even  during  large  storms.  This  is  somewhat  surprising,  considering  the 
results  of  O’Brien  and  McPherron  [2000]  who  conclude  that  the  Dst  activity  can  be 
almost  completely  predicted  from  the  solar  wind  driver  parameters. 

The  correlation  of  the  ion  fluxes  with  the  F10.7  fluxes  is  weaker  than  one  would 
think,  when  many  previous  results  have  shown  that  the  ionospheric  ion  outflow  is 
well-correlated  with  the  ion  fluxes  [e.g.,  Young  ct  a/.,  1982:  Yau  at  al...  1988:  Moor <  et. 
al .,  1999].  Figure  10a  shows  four  frames  taken  from  the  Yau  ct  al.  paper,  which  give  the 
average  response  of  the  ionospheric  ion  outflow  as  measured  by  tin1  Dynamics  Explorer 
(DE)  1  spacecraft  during  1981-  1986.  The  figure  shows  that  the  H+  fluxes  increase  as  a 
function  of  activity  level  (Kp  or  Dst).  but  there  is  no  change  for  different  F10.7  flux 
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Figure  10a.  Ionospheric  ion  outflow  rate  of  H+  and  0+  ions  at  0.01-17  keV  as  a  function  of  the 
magnetic  Kp  index  and  Dst  index  for  different  ranges  of  FI  0.7  values  as  measured  by 
Dynamics  Explorer  1  during  1981-1986  (from  Yau  etal.  [1988]). 


levels.  On  the  other  hand,  the  0+  fluxes  show  clear  differences  when  the  data  are 
separated  according  to  the  F10.7  flux:  separate  linear  fits  can  be  made  which  indicate 
flux  levels  almost  an  order  of  magnitude  higher  during  high  F10.7  flux  than  during  low 
solar  activity. 

Figure  10b  shows  the  average  energetic  H+  and  0+  fluxes  measured  in  the 
magnetosphere  at  similar  intervals  of  F10.7.  A  similar  trend  is  evident:  H+  fluxes  show 
little  correlation  with  F10.7  fluxes,  but  the  0+  signature  is  much  less  clear.  First,  the 
trends  for  increasing  activity  are  not  as  linear  as  they  arc  for  the  Yau  et  al.  result, 
and  the  differences  between  the  high  solar  activity  and  low  solar  activity  are  not  as 
clear.  Especially,  for  high  activity  levels,  the  values  are  not  well-separated  from  each 
other.  This  indicates  that  there  are  other,  more  complicated  processes  which  control  the 
acceleration  of  the  ions  to  the  ring  current  energies.  This  would  suggest  that  even  for  a 
perfect  data  set,  the  correlation  analysis  might  not  lead  to  better  results.  It  appears  that 
some  significant  effects  are  introduced  in  going  from  the  polar  outflow  region  through 
the  equatorial  plasma  sheet  and  into  the  ring  current  proper. 

Observations  of  lower-energy  ions  at  geosynchronous  orbit  [  Young  et  al.,  1982: 
Stockholm  et  al.,  1989]  show  much  better  correlations  between  ion  densities  and  the 
F10.7  fluxes  than  those  obtained  in  this  study.  The  observations  made  by  the  European 
Space  Agency’s  Geosynchronous  spacecraft  (GEOS  1  and  GEOS  2)  Ion  Composition 
Experiment  used  in  those  studies  covered  an  energy  range  of  1  -  15  keV,  i.e..  energies 
typical  for  particles  that  have  been  accelerated  from  the  auroral  region  ionosphere  ro 
the  magnetotail  near  geosynchronous  orbit.  On  the  other  hand,  the  POLAR/MICS 
measurements  used  in  this  paper  that  extend  to  much  higher  energies  and  to  lower 
L  shells.  These  measurements  are  in  an  energy  regime  where  additional  acceleration 
processes  have  occurred  after  the  initial  acceleration  from  the  ionosphere,  and  in  a 
spatial  region  where  they  have  been  transported  from  the  more  distant  magnetotail. 

The  ring  current  ions  enter  the  inner  magnetosphere  mainly  by  radial  diffusion 
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Figure  10b.  Ring  current  H+  AND  0+  ion  fluxes  at  1-200  keV  as  a  function  of  the  magnetic  Kp  index 
and  Dst  index  for  different  ranges  of  FI  0.7  values  as  measured  by  POLAR/CAMMICE/ 
MICS  during  1996-1999. 


from  the  tail.  Kremser  et  al.  [1993]  used  data  from  the  CHarge-Energy-Mass  (CHEM) 
instrument  onboard  the  Active  Magnetospheric  Particle  Tracer  Explorer  /  Charge 
Composition  Explorer  (AMPTE/CCE)  satellite  together  with  a  model  for  radial 
diffusion.  They  obtained  a  good  equilibrium  between  radial  diffusion  and  charge 
exchange  losses  for  helium  inside  of  about  7  Re-  On  the  other  hand,  Sheldon  and 
Hamilton  [1993]  examined  only  the  very  quiet-time  ring  current,  and  concluded  that  a 
pure  radial  diffusion  model  is  not  able  to  account  for  the  observed  high  fluxes  at  L  <  4. 
They  assumed  that  this  discrepancy  was  due  to  lack  of  particle  scattering  caused  by 
electric  field  fluctuations  in  their  model.  The  complex  processes  that  occur  during  the 
inward  transport  of  the  ions  are  a  likely  cause  to  the  large  scatter  in  the  correlations 
obtained  in  this  study. 

According  to  the  Dessler-Parker-Sckopke  relation  [em  Dessler  and  Parker,  1959; 
Sckopke,  1966],  the  energy  of  the  ring  current  particles  is  proportional  to  the  ground 
disturbance  of  the  geomagnetic  field  with  a  proportionality 


A B  __  2W 


(3) 


where  A B  is  the  field  decrease  due  to  the  ring  current,  B0  is  the  average  equatorial 
surface  geomagnetic  field,  W  =  f  wdV  is  the  total  energy  of  the  ring  current  particles 
or  the  energy  density  integrated  over  the  ring  current  volume,  and  Wm  ~  1018  J 
is  the  energy  in  the  Earth’s  dipole  field  above  the  surface.  Assuming  that  the  ring 
current  volume  is  approximately  constant,  the  ground  magnetic  disturbance  should  be 
proportional  to  the  energy  densities  measured  by  POLAR/MICS.  Figure  11  shows  the 
correlation  between  the  Dst  index  and  the  summed  energy  densities.  The  correlation 
coefficient  is  about  the  same  as  it  is  for  the  individual  species.  Thus,  different 
partitioning  of  the  energy  density  between  the  different  ion  species  is  not  responsible  for 
the  scatter  in  the  correlations  shown  in  this  paper. 

To  examine  the  scatter  produced  by  the  fact  that  not  all  measurements  were  made 
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near  the  peak  activity,  the  storm  data  set  was  limited  to  those  where  the  measurement 
was  made  within  2  hours  of  the  minimum  of  Dst.  This  limited  the  data  set  further 
to  only  12  events,  but  the  correlation  analyses  yielded  practically  identical  results. 
Assumption  of  the  constant  volume  is  another  source  for  error;  however,  inclusion  of 
a  variable  volume  would  have  to  be  made  on  an  event-by-event  basis,  and  is  left  as  a 
future  study. 

Part  of  the  scatter  in  the  correlations  of  course  comes  from  the  definition  of  Dst 
itself:  instead  of  being  a  pure  measure  of  the  ring  current,  it  responds  to  all  other 
magnetospheric  current  systems  as  well:  The  magnetopause  current  is  often  corrected 
by  the  solar  wind  dynamic  pressure  effects  [Burton  et  al,  1975],  the  tail  current  has 
been  shown  to  contribute  about  25  nT  close  to  the  time  of  the  substorm  onset  when 
the  cross-tail  current  is  strong  [Turner  et  al.,  2000a].  Furthermore,  the  ring  current 
is  strongly  asymmetric,  which  then  leads  to  different  flux  values  depending  on  the 
longitude  the  measurements  were  made.  Lastly,  the  effects  caused  by  induced  currents 
in  the  ground  add  about  30%  to  the  measured  variation  [Turner  et  al.,  2000b].  Given 
all  these  factors  which  affect  the  Dst  measurements  in  a  way  different  for  each  event, 
one  would  not  expect  to  have  a  perfect  correlation  between  the  Dst  and  the  ring  current 
ion  fluxes. 

It  is  interesting  to  note  that  the  0+  fluxes  give  best  correlations  with  the  peak  of 
Dst  rather  than  the  specific  value  of  the  Dst  at  the  time  of  the  flux  measurements.  This 
would  indicate  that  strong  enhancement  of  0+  requires  strong  substorm  activity  during 
the  main  phase,  which  could  increase  the  outflow  from  the  ionosphere.  However,  even 
the  correlation  with  peak  substorm  activity  (minimum  of  AL  index)  did  not  show  very 
high  correlation.  A  more  thorough  evaluation  of  the  substorm  activity  to  the  outflow  is 
left  as  a  future  study. 

To  conclude,  in  this  paper  we  demonstrate  clear  changes  in  the  ring  current 
composition  and  intensity  as  the  solar  conditions  changed  from  the  minimum  toward 
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maximum  activity  during  the  three  years  of  POLAR  operations.  POLAR  was  launced 
in  a  very  fortuituous  time  to  record  the  low  fluxes  and  quiet  magnetosphere  near  the 
solar  minimum,  and  has  continued  to  produce  data  until  the  solar  maximum.  From  the 
trends  here,  one  can  conclude  that  there  are  distinct  effects  that  are  associated  with  the 
long-term  solar  variability,  which  is  indicative  of  a  solar  cycle  dependence  of  the  inner 
magnetosphere  trapped  population. 

The  correlation  analyses  shown  here  demonstrate  that  the  processes  governing 
the  ring  current  intensity  are  complex  functions  of  the  solar  wind  driver  as  well  as 
of  the  intramagnetospheric  processes  such  as  substorms.  Therefore,  the  next  step  in 
understanding  the  processes  will  have  to  involve  detailed  analyses  of  each  of  the  storm 
events  separately,  addressing  the  temporal  sequence  for  each  case  separately.  This  is  left 
as  a  future  study. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  special¬ 
izing  in  advanced  military  space  systems.  The  Corporation's  Laboratory  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staff’s  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analy¬ 
sis,  solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid 
state  laser  design,  micro-optics,  optical  communications,  and  fiber  optic  sensors;  atomic 
frequency  standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation 
and  beam  control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evalua¬ 
tion,  battery  electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  fail¬ 
ure  analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  space  environment 
effects  on  materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and 
structural  control;  lubrication  and  surface  phenomena. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing 
using  atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
infrared  surveillance,  imaging,  remote  sensing,  and  hyperspectral  imaging;  effects  of  solar 
activity,  magnetic  storms  and  nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and 
magnetosphere;  effects  of  electromagnetic  and  particulate  radiations  on  space  systems; 
space  instrumentation,  design  fabrication  and  test;  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific  chemical 
reactions  and  radiative  signatures  of  missile  plumes. 

Center  for  Microtechnology:  Microelectromechanical  systems  (MEMS)  for  space 
applications;  assessment  of  microtechnology  space  applications;  laser  micromachining; 
laser-surface  physical  and  chemical  interactions;  micropropulsion;  micro-  and  nanosatel¬ 
lite  mission  analysis;  intelligent  microinstruments  for  monitoring  space  and  launch  sys¬ 
tem  environments. 


Office  of  Spectral  Applications:  Multispectral  and  hyperspectral  sensor  development; 
data  analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral 
imagery  to  defense,  civil  space,  commercial,  and  environmental  missions. 


